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On the composition of finite rotations in 


Alex Goldvard* and Lavi Karpl 


Abstract 

We achieve compositions rules for the geometric parameters of the 
composed rotations, which is in a certain sense analogous to the well 
known Rodrigues formula. We also obtain a necessary and sufficient 
condition for a composition of two simple rotations in to be a simple 
rotation. 


1 Introduction 

In this note we are studying compositions of finite rotations in the Euclidean 
space the four dimensional real vector space with the standard scalar 
product. 

There are several significant differences between rotations in and in the 
three dimensional space. The group of rotations in E^ essentially comprises 
one type, that is, a rotation about an axis, while in E^ there are three 
types of rotations: (i) A simple rotation, it leaves a plane (two dimensional 
subspace of E^) point-wise fixed and induces a two dimensional rotation in 
the orthogonal plane; (ii) A Clifford translation (an isoclinic rotation), here 
each vector in E^ turns through the same angle; (hi) A double rotation, here 
E^ is decomposed into two orthogonal planes and points in the first plane 
rotate through an angle a, while points in the second plane rotate through 
an angle, f a. 

In this note we derive two properties of the compositions of rotations in 
E^. The first one gives analytical and geometrical characterizations of the 
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subgroup of simple rotations. The second result deals with compositions 
formulas for double rotations, that is, we obtained formulas that enable the 
calculations of the orthogonal planes and angles of the composed rotation 
in terms of the corresponding characteristics of each one of the rotations in 

Our main tool is the quaternionic representation 

X 1-^ axb (1) 

of rotations in E^. Cayley obtained formula (1) in 1855 by means of labo¬ 
rious calculations [2]. Many years later, in 1945, Coxeter [4] gave a very 
elegant proof of the above formula by employing compositions of reflections 
to represent rotations. He also derived the geometric characterization of the 
rotation ( 1 ), that is, he computed the orthogonal planes and angles via the 
quaternions a and b. But there are several cases in which Coxter’s formulas 
are not valid. We shall fill the gap in this note and complete the geometric 
characterization in those cases. 

The outline of the paper is as follows. In the next section we shall give a brief 
acconnt of Coxeter’s approach and treat the above mentioned special cases. 
Section 3 deals with the composition of simple rotations. In section 4 we 
calculate planes and angles of a composed rotation and exhibit an example. 
We end with section 5 where we discuss the results. 


2 Representation of rotations and reflections by 
quaternions 

We first fix the notations. A quaternion x is expressed as x = XQ+Xii+X 2 j + 
x^k = Sx + Vx, where Sx = xq is called the real part, Vx = xii + X2j + x^k 
is the vector part, and the coefficients xo,xi,X 2 and 0:4 are real numbers. 
The quaternions multiplication is defined by the table 

= k^ = — 1 , ij = —ji = k,jk = —kj = i,ki = —ik = j. 

We denote by El the set of all quaternions. Since a quaternion is completely 
defined by its four real components, it determines a point (vector) in the four 
dimensional real space Similarly, Vx = Xii + X 2 j Vx^k can be identified 
with a point (vector) in three dimensional space The multiplication 
formula 


xy = SxSy — Vx ■ Vy + SxVy + SyVx + Vx x Vy, (2) 
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is a consequence of the distributive property of the quanternionic algebra. 
Here 

Vx ■ Vy and Vx x Vy are the scalar cross products in respectively. 

Let X be a quaternion, then x = Sx — Vx is the conjugate quaternion, and 
Nx = XX = XX = XQ + xf+ X 2 + X 3 is the square of the norm. When Nx = 1, 
then it is called a unit quaternion. Any unit quaternion a can be written in a 
polar form a = cos a+p sin a, where p = Pii+P 2 j+P 3 k and satisfies = — 1 . 
The vector p is called a pure unit quaternion. A subspace which is spanned 
by two quaternions a and b is denoted by Sp{a, b} = {S 10 + S 26 : si, S 2 € M}. 

The advantage of quaternionic representation relies on two simple properties 
of quaternions. The function /(x) = ox is linear and N{ax) = NaNx. 
Therefore, if a is a unit quaternion, then the multiplication ax represents 
an orthogonal transformation in E^. Moreover, for any unit quaternion x, 
(ox) • X = S'(axx) = S{a), hence the cosine of the angle between ax and x is 
equal to S{a) for any x £ E'^. Thus any vector in E^ rotates with angle a, 
where cos a = S{a). Such type of rotations does not exist in E^, and it is 
called a left Clifford translation. It is evident that the set of all left Clifford 
translations forms a subgroup of the rotations’ group in E^. Similarly, one 
may consider a right Clifford transformation x i—>• xb, where 6 is a unit 
quaternion. Now, the composition of left and right Clifford translations 
results in the formula (1). As we already mentioned, Coxeter proved that 
this formula comprises all rotations in E^ by using reflections with respect to 
hyperplanes. For readers’ convenient we recall his description of reflections 
and rotations. 

Reflections: 

Notice that any two quaternions y and x are orthogonal in E^ if and only if 
xy + yx = 0. Hence if Ny = 1, then the transformation 

X —yxy 

is a reflection with respect to a hyperplane having y as its normal (see [4, 
§5]). In what follows we denote this reflection by Ry. 
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Simple rotations: 


A rotation in the plane through the angle a and all the points in the or¬ 
thogonal plane remain fixed. Its representation is given by 

OL 

X !->■ axb, Na = Nb = 1, Sa = Sb = cos —. (3) 

It is shown in [4] that there are two unit quaternions y and z such that a = 
zy, b = yz, and the simple rotation (3) is the composition of the reflections 
Ry{x) = —yxy and Rz{x) = —zxz, that is, 

X 1-^ axb = {zy)x{yz). (4) 

Obviously, x >->■ {zy)x{yz) = x, whenever x is orthogonal both to y and 
z. Hence, the point-wise invariant plane of the simple rotation (4) is the 
intersection of the hyperplanes 

3 3 

Hi := {x G M; ^ = 0}, II 2 := {x G H; ^ z^Xf, = 0}. 

y,=0 y,=0 


Clifford translations: 

Each vector in turns through the same angle. Let a be a unit quaternion, 
then X e->■ ax is a left Clifford translation and x 1 —)• xa is a right Clifford 
translation. The angle a of both types satisfies the equation cos a = Sa. 

Writing a in a polar form a = cos a-tpsina, and consifer the left translation 
X ax, then 


X !->■ xcosa+ pxsina, px >-)• pxcosa — xsina. (5) 

Hence the plane Sp{x,px} is an invariant subspace of this transformation, 
and consequently Clifford translations have infinite number of invariant two 
dimensional subspaces. Similarly, for the right Clifford translation x xb, 
b = cos/3 + (/sin/3, the planes Sp{x,x(/} are invariant subspaces. These 
subspaces are completely defined by the pure quaternions p and q of the 
Clifford translations. In addition, we observe the similarity between (5) and 
the formula for rotations in E^. 
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Double rotations: 


Two orthogonal planes rotate simultaneously with two different angles. Let 
p and q be two pure quaternions such that p^ = = —1, then a double 

rotation has the representation 

X !->■ (cosa +psina)x(cos/3 + qsm/3) =: f{x). ( 6 ) 

The rotations’ angles are a±/3. In case that cos a = cos/3, then the formula 
( 6 ) represents a simple rotation. 

Note that if cos a, cos/3 / 0, then we can set p = ptana, q = gtan/3 and 
rewrite formula ( 6 ) as 

f{x) = cosa{l+p)x{l + q)cos/3. (7) 

We will use this representation in Section 4. The quaternions p and q are 
called Gibbs vectors. We shall call (7) Gibbs representation of the rotation. 

The orthogonal planes of the double and simple rotations: 

Goxeter proved that the orthogonal planes of the rotation ( 6 ) are 

Hi = Sp{p - 1+PQ'} and U 2 = Sp{p + q,l - pq} . 

Points in the plane Hi rotate through the angle a + /3 and in the plane 112 
through the angle a — 13, see [4, §9]. It follows that if /3 = a, then ( 6 ) is 
a simple rotation with the points-wise invariant plane 112 , while Hi is the 
fixed-plane when f3 = —a. 

If one of the quaternions p =F g, 1 ± g is zero, then Coxeter’s characterization 
of the invariant planes is not longer valid and we have to treat these cases 
separately. First of all note that the equalities p = ±g and pq = =pl are 
equivalent and, in fact, there are only two special cases: q = ±p. It is 
readily follows from formula ( 6 ) that in these cases 

/(l)GSp{l,p} and /(p) G Sp {l,p} . 

Therefore Ai = Sp{l,p} is one of the rotating planes. Since the second one 
A 2 , is the orthogonal complement of Ai, we conclude that 

A 2 = {x G E[;p • X = = 0} 

The pointwise invariant planes of simple rotations when q = ±p can be found 
by the following considerations. q = p, then by /(I) 7 ^ 1 , and therefore 
A 2 is the fixed plane. While if g = —p, then /(I) = 1, and hence Ai is the 
fixed plane. 
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3 Composition of simple rotations in 


Let / ; X axb and g : x cxd be two simple rotations and consider the 
composition go f ■. x ^ {ac)x{bd). We ask under which conditions go f will 
be again a simple rotation. 

It follows from (3) that the rotation g o f is simple if and only if S{ca) = 
S{bd). According to the multiplication formula ( 2 ), this condition is equiv¬ 
alent to 

ScSa -Vc-Va = SbSd - Vb ■ Vc. 

Since both / and g are simple rotations, S(a) = S{b) and 5(c) = 5(d), hence 
the above condition results in 


Vc-Va-Vb-Vd = 0. ( 8 ) 

We now recall that each rotation is decomposed into two reflections by for¬ 
mula (4). Hence, there exists four unit quaternions y,z,u,w such that 

a = zy, b = yz, c = wu, d = uw. 


Hence 

Vc-Va = V{wu)-V{zy) 

= [—wqVu -|- uqVw + Vu X Vw] ■ [—zoVy + yoVz + Vy x Vz] 

and 

Vb-Vd = V{yz) ■ V{uw) 

= [yoVz — zqV y -VyxVz]- [uqVw — wqVu — Vu x Vw] . 

Thus (8) becomes 

Vc-Va-Vb-Vd = -2 {yoVz • {Vu x Vw) - zoVy • {Vu x Vw) 

+uoVw ■ {Vy X Vz) — wqVu ■ {Vy x Vz)} (9) 

= —2 det[y, z, u, w] = 0, 

where [y, z, u, rc] denote the matrix with the column vectors y, z, u and w. 
These vectors are the normals of the hyperplanes of the reflections Ry,Rz, Ru 
and Rw that decompose the rotations. Thus we have obtained the following 
theorem. 
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Theorem 1. Let f : x axb and g : x cxd be two simple rotations 
such that f is the composition of the reflections Ry with R^, and g is the 
composition of the reflections Ru with Ryj. Then the composition g o f is a 
simple rotation if and only if the four normals of the reflections are linearly 
dependent. 

The point-wise invariant planes of the rotations / and g are 

3 3 

11 / := {x £ M-,'^yyXy = '^ZyXf, = 0 } 

/^=0 fj .=0 


Ug := {x G H; ^ UgXg = ^ WgXg = O} . 

/^=0 ^ 1=0 

Obvionsly, if 11/ n 11^ {0}, then the composition go f is simple. However, 

the converse implication is not trivial. It follows from the above Theorem, 
that if the composition is simple, then the four normals are linearly depen¬ 
dent, which implies that the intersection H/riHg contains a non-zero vector. 
So a consequence of the Theorem 1 is: 

Corollary 1. Let f : x axb and g : x cxd be two simple rotations, 
and H/ and H^ be the point-wise invariant planes of the rotations f and g 
respectively. Then the composition h = go f is a simple rotation if and only 

^/^/n^3/{o}. 

Remark 1. When the paper was written the authors were not aware that 
in 1890 Cole proved Corollary 1. Cole essentially worked out Cayley’s pa¬ 
rameters of the subgroup and this method demands laborious computations. 
He also obtained formulas for the orthogonal planes and the angle of the 
composed rotation, but they are complicated and lengthy, and inapplicable 
even in some very simple examples. 


4 The angles and planes of the composed rotation 

In this section we consider the following problem: Does there exist formu¬ 
las that link the orthogonal planes and angles of two rotations to those of 
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composed rotation. We start with double and simple rotations. So suppose 
that 


f{x) = axb = cosai(l + pi)x{l + qi) cos /3i ( 10 ) 

g{x) = cxd = cosa 2 (l + P 2 )x{l + 92 ) cos /32 (H) 

are two double rotations such that 

cos «!, cos/3i 7 ^ 0 and cos 02 , cos /32 / 0 . ( 12 ) 

Recall that right hand side of (10) and (11) are Gibbs representation of the 
rotations / and g. Consider the composition / followed by g 

h{x) := {g o f){x) = {ca)x{bd) = cosa(l + p)x{l + q) cos/3. 

Then in accordance to the multiplication formula (2), 

cos a = S{ca) = cos ai cos 02(1 — Pi ■ P 2 ) (13) 


and 
Assuming 
we obtain that 


p sin a = V(ca) = cos ai cos 02 {P 2 + Pi + P 2 Pi)- 
cos a / 0 , 


p tan a = p = 


P 2 +P 1 +P 2 X Pi 


^-P2-Pi 

In a similar manner, 

cos P = S{bd) = cos /3i cos / 32(1 — qi ■ ^ 2 ), 
qsin/3 = V{bd) = cos/3i cos/ 32 (gi + q 2 + Qi x ^ 2 ), 


and if 

then 


5 tan (3 = q = 


cos /3 7 ^ 0, 

qi + q2 + qix q2 


(14) 

(15) 

(16) 

(17) 

(18) 


1-qi -^2 

Thus in a view of those computations, we have obtained the following the 


orem. 




Theorem 2. If the conditions (12), (If) and (17) hold, then the compo¬ 
sition of the rotations ( 10 ) with ( 11 ), h := g o f, is a rotation through the 
angles ah /3 about the planes spanned by quaternions ptana ± qtan(3 and 
tana tan/3=FM; where a, fi, p and q are given by (13) and (16), (15) and 
(18) respectively. 


Note that the angels of the composed rotation are independent of the order 
of the composition, that is, g o f and fog result in the same angles. While 
the cross-product in (15) and (18) depend upon the order. 

Remark 2. The expressions (15) and (18) resemble the known Rodrigues 
formula for the composition of rotations in (see e.g. [7, p. 58, 71]). 
That is, if f is a rotation through an angle ai about an axis ci, and g is a 
rotation through an angle about an axis C 2 , then their composition g o f 
(f followed by g) is a rotation through an angle a about an axis c, where 


C2 + Cl + C2 X Cl 
1 - C2 • Cl 


(19) 


c := ctan^jc'i := ci tan ^ and C 2 := C 2 tan^. The similarity of the for¬ 
mulas (15), (18) and (19) is not accidental. A rotation in E^ about an axis 
c through an angle a has the quanternionic representation x i—)• axa, where 
a = cos § + csin Therefore, the formula (19) is obtained in the same way 
as we do it for the composition of two double rotations in E^. The vector c 
is called Gibbs vector after J. W. Gibbs who invented it in 1881 in order to 
achieve formula (19). In fact this formula was obtained by Olinde Rodrigues 
in I 84 O, see [6] for historical survey. 

Remark 3. Note that for a rotation in E^, c = ctan^. Hence the dis¬ 
continuity occurs when a = Htt, which means that the composed rotation is 
a reflection. For rotations in E^ there are many possibilities for disconti¬ 
nuity. For example, in the case of a composition of simple rotations, where 
cos a = Q in (13), then it results in a reflection in one plane, and the identity 
in the orthogonal plan. 


Remark 4. To the authors’ knowledge Fedorov was the first one who dis¬ 
covered compositions rules of four dimensional rotations in 1958. Fedorov 
interest was in applications of rotations to relativity and quantum physies, 
hence he considered Lorentz transformations in Minkowski space. He ob¬ 
tained analogous formulas to (13) and (16), (15) and (18) in the Minkowski 
space and E'^. Fedorov showed that the group of rotations in E^ can be 
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parametrized by two arbitrary three dimensional vectors, and by writing ex¬ 
plicitly the matrices of rotations in via those vectors he derived the com¬ 
positions formulas ([5], ^24)- So Fedorov basically used matrices approach 
as opposed to the quaternions algebra. Fie did not consider the geomet¬ 
rical characterization (orthogonal planes and angles). Recently Fedorov’s 
ideas on using Gibbs parametrization of a rotation (it is also called a vec¬ 
tor parametrization) became quite widespread, see for example [1],[8] for 
further discussions. 


We turn now to the composition of two Clifford translations of the same type. 
So suppose f : X ax and g : x bx are two left Clifford translations, and 
consider the composition / followed hy g: h = go f ■. x cx, where c = ba. 
Writing the quaternions by means of Gibbs vectors, that is, a = cosai(l + 
Pi), b = cosa 2 (l + P 2 ) and c = cosa(l + p), then by the quanternionic 
algebra we get that 

. Pl+P2+P2^ Pi 

P= -- 

I - Pi ■P2 


and 


cos a = cos cos a2(l — Pi • P2)- 


Thus in view of formulas (13), (15), (16), (18), (19) and the two last ones, 
we may say that in a certain sense all type of rotations obey the same 
composition rule. 


4.1 Example 


The following example illustrates the consequences of Corollary 1 and The¬ 
orem 2. Let 


1 + il + j 1 -\- j 1 -\- k 

t : X — ^x — ^ and q : x — ^x—^ 

V2 V2 V2 V2 


be two simple rotations. Then 11^ = Sp {z — j, 1-|-A:} is the fixed-points 
plane of the rotation / and lig = Sp {j — k,l -\- {} of g. We see that n 
llg 7 ^ {0}, hence by Corollary 1 the composition g o f should be a simple 
rotation. Simple calculations show that 


h 


go f -. xe^ 


1 + ^ + j 
2 


k 1-\-i-\- j -\-k 

—X - - -=: axb. 


Since Sa = Sb, the composition gofis simple. We shall now verify formulas 
(13), (15) and (18). In order to it we represent these rotations by the 
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Gibbs vectors as in (7). So let Pf,qf,Pg,qg,Ph,Qh the Gibbs vectors that 
corresponding those rotations, and let af, ag, be angles of Gibbs vectors 
(p = plan a, where p is a pure unit quaternion). Then 


/ : X 1 -^ (^cos ^ (1 + *)) X (^{l + j) cos , 

g : (cos ^ (1 + j)) x (^(1 + k) cos 


and 


h \ X ^ ^cos — (1 + i + j — k^ X ^(1 + i + j + k) cos —^ . 


By (13), 


cos ah = cos a/- cos ag {1 — {i ■ j)) = (cos —)^ = 


( 20 ) 


Thus ah = ^, which obviously coincides with (20). Now, in accordance with 
formulas (15) and (18), 


Ph 


j + i + j xi 
I-j - i 


i + j — k and qh 


j + k + j X k 

1- j -k 


i + j + k. 


5 Conclusions 

There are various methods to represent rotations in lE^. In this article we 
are using the quaternionic algebra approach. It enables us to obtain the 
characterization of a nontrivial subgroup of the simple rotations. Another 
advantage is the quaternionic formula (6) for double rotations that describes 
in a very intrinsic way the geometric parameters. Relying upon that formula 
we had achieved the compositions rules by means of very simple calculations. 
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